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ABSTRACT 


Analysis  of  the  procedures  recommended  by  Marshall,  Springer 
and  Rodean  (1979)  for  correction  of  explosion  magnitude  due  to  pP  and 
absorption  shows  that  the  techniques  are  needed  and  are  generally  satisfactory. 
However,  the  pP  correction  could  be  made  more  accurate  and  operationally 
applicable,  and  the  absorption  correction  needs  further  checking  and  perhaps 
a  different  approach.  A  new  attack  on  the  general  problem  is  recommended  in 
order  to  properly  treat  problems  of  focus ing-defocusing,  crustal  amplification, 
clipping  and  noise  data,  and  pP  and  absorption  as  mentioned  above. 
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WWSSN  waveforms  and  short-period  measurements  of 
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seconds.  Full  and  partial  pP  reflection  coefficients. 


3- 


LIST  OF  FIGURES 


Figure  No. 


Title 


Page 


1  Observed  relationship  between  Pn  and  1°Riq(A/T)  residuals,  9 
comparison  with  calculated  Amb  corrections  and  0a  for 

T  =  0.75s.  (from  Marshall  et  al,  1979). 

2  LRSM  magnitude  residuals  from  Booth  et  al  (1974)  with  10 

crustal  corrections  due  to  Der  et  al  (1979)  as  a  function 

of  Pn  velocity  as  from  Marshall  et  al  (1979).  The  curve 
drawn  is  from  Marshall  et  al  adjusted  down  by  0.12  to 
allow  for  the  crustal  corrections. 

3  Magnitude  residuals  of  North  (1977)  plotted  as  a  function  12 

of  Pn  velocitv  due  to  Marshall  et  al  (197°).  The  super¬ 
imposed  line  is,  again,  due  to  Marshall  et  al  (1979). 

4  Mean  biases  for  continental  US  stations  versus  13 

(a)  sub-station  Pn  velocity  from  Herrin  (19f2); 

(b)  station  P  travel-time  anomalies  from  Hales  and 
Herrin  (1972). 

5a  Variation  of  the  pP-P  interference  effect  with  explosion  16 

depth  (from  Marshall,  Springer  and  Rodean,  1979); 

5b  The  magnitude-correction  term  DC(T)  given  as  a  function  16 

of  the  depth-correction  ratio  DCR.  The  DC(T)  is  added  to 
magnitude  to  give  the  final  magnitude  mQ  (from  Marshall, 
Springer  and  Rodean,  1979); 

5c  Magnitude  correction  as  a  function  of  DCR  =  (pP  delay/  17 

c  phase  period)  from  waveforms  in  Appendix  I.  Perfect 
pP  reflection  (apP  *>  -1.0),  WWSSN  SP  response, 
von  Seggern-Blandford  granite  reduced  displacement 
potential.  Correction  is  to  be  added  to  achieve  maximum 
mb  possible  with  perfect  pR  enhancement  for  each  yield; 

5d  Magnitude  correction  as  a  function  of  DCR  ■  (pP  delay/  17 

c  phase  period  from  waveforms  in  Appendix  I.  Variable 
pP  reflection  (-0.5(1  +  exp(Tf2)),  WWSSN  SP  response, 
von  Seggern-Blandford  granite  reduced  displacement 
potential.  Correction  is  to  be  added  to  achieve  maximum 
mb  possible  with  perfect  enhancement  for  each  yield. 

6  RKON  station:  Residuals  of  log  (a)  with  respect  to  a  c  18 

as  a  function  of  scaled  2-way  pP  travel  time.  The  solid 
line  is  hand-drawn  as  a  subjective  estimate  of  the  main 
data  trends  (from  Blandford,  1976b). 


-4- 


Figure 

7 

8 


9 


10 


11 


12 

13 


LIST  OF  FIGURES  (Continued) 

No.  Title  Page 

Residuals  of  log  (a)  at  HNME  with  respect  to  c  as  a  19 

function  of  scaled  2-way  pP  travel  time.  The  solid 
line  is  hand-drawn  as  a  subjective  estimate  of  the 
mean  data  trends  (from  Blandford,  1976b). 

Theoretical  amplitude-yield  curves  for  t*  =  0,  0.1,  20 

0.2,  0.4,  0.6;  granite,  one-half  maximum  peak-to- 
trough  amplitude  for  signal,  corrected  for  LRSM  instru¬ 
ment  response  at  measured  period,  T,  and  divided  by 
T,  with  surface  reflection.  pP  delay  equal  to  0.12Y^/3 
sec,  with  Y  in  kt.  From  Blandford  (1976b). 

Theoretical  amplitude : yield  for  t*  =  0.4,  granite;  22 

c  phase  corrected  for  WWSSN  response  at  measured 
period  T,  and  divided  by  T,  with  surface  reflection 
coefficient  apP(f)  =  -0.5(1  +  exp(-f^)),  pp  delay 
equal  to  .05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y  =  75, 

150,  300  kt.  From  waveforms  in  Appendix  I. 

Theoretical  amplitude: yield  for  t*  =  0.4  granite;  c  phase  23 
corrected  for  WWSSN  response  at  measured  period  T,  and 
divided  by  T,  with  surface  reflection  coefficient  apP(f)  = 

-1.0  pP  delay  equal  to  .05,  .1,  .15,  .2,  .3, 

.4,  .6  sec  for  Y  =  75,  150,  300  kt.  From  waveforms  in 
Appendix  I. 

Theoretical  curves  of  magnitude  and  period  as  a  function  25 

of  t*  for  several  definitions  of  magnitude.  Granite,  WWSSN 
response  surface  reflection  coefficients  apP(f)  » 

-0.5(1  +  exp  (-f^)),  pP  delay  equal  to  0.15  sec,  Y  =>  150  kt. 

C  is  c  corrected  for  instrument  response  at  period  T,  c/g, 
afid  m  ub  c/GT.  From  waveforms  in  Appendix  I. 

.26 

_  Theoretical  amplitudervield  curves  for  t*  ■  0.4  granite, 
phase  with  surface  reflection  coefficient  apP(f)  *  -1.0, 
pP  delay  equal  to  .05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y 
■  75,  150,  300  kt.  From  waveforms  in  Appendix  I. 

Theoretical  amplitude:yield  curves  for  t*  *  0.4  granite,  27 

£  phase  with  surface  reflection  coefficient  apP(f)  -  -1.0, 
pP  delay  equal  to  .05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y 
*  75,  150,  300  kt.  From  waveforms  in  Appendix  I. 


-5- 


INTRODUCTION 


A  recent  paper  by  Marshall,  Springer  and  Rodean  (1979),  hereafter 
referred  to  as  MSR,  has  presented  a  comprehensive  procedure  for  calculating 
yields  from  short-period  body  wave  magnitude.  In  that  report,  they  give 
procedures  for  correcting  the  magnitude  for  effects  of  source  and  receiver 
absorption,  pP,  and  station  effects.  The  principal  objective  of  this 
report  is  to  evaluate  their  procedures  for  absorption  and  pP  corrections. 


Attenuation  Correction 

The  basic  approach  in  MSR  is  to  determine  mantle  Q  under  inaccessible 
locations.  While  the  magnitude  and  spectral  measurements  need  for  deter¬ 
mining  0  are  not  available  in  most  regions,  detailed  crustal  surveys  that 
include  values  of  P^  velocities  have  been  published  for  areas  of  the  Soviet 
Union,  China,  South  Africa,  etc.  MSR  use  the  apparently  good  correlation 

between  P  velocities  and  the  residuals  of  Booth  et  al.  (1974)  to  predict 
n 

m^  residuals  for  the  rest  of  the  world.  While  the  idea  is  attractive,  more 
work  is  needed  to  establish  the  validity  of  such  a  correlation  worldwide.  The 
residuals  of  Booth  et  al.  (1974)  are  not  corrected  for  crustal  effects.  It 
has  been  found  that  such  corrections,  while  somewhat  inaccurate,  can  be  used 
to  reduce  the  variance  in  station  biases  (Der  et  al.,  197**).  Beside  the 
lack  of  crustal  corrections,  there  is  considerable  variation  among  the  various 
studies  of  magnitude  biases  for  the  same  set  of  stations  and  the  neat  correla¬ 
tion  that  is  displayed  by  MSR  may 'not  exist.  Another  possible  objection  to 
MSR's  procedure  is  that  the  physical  basis  for  causal  interrelationship  between 
Pn  velocity  and  mantle  attenuation  is  somewhat  shaky.  The  Pn  velocity  is  deter¬ 
mined  by  the  uppermost  mantle,  the  so  called  "lid  zone",  while  the  magnitude 
bias  is  caused  by  losses  in  the  total  thickness  of  the  upper  mantle.  Never¬ 
theless,  there  may  be  a  causal  link;  portions  of  the  mantle  with  low  0  are 
probably  hot,  which  causes  lowering  of  the  P^  velocity.  However,  composi¬ 
tional  changes  are  also  necessary  to  explain  the  range  of  P^  velocities  under 
the  United  States. 

In  order  to  test  the  validity  of  0  corrections  derived  from  P  data  we  have 

n 

corrected  the  residuals  of  Booth  et  al  for  crustal  effects  using  the  results  in 
Der  et  al.  (1979)  and  plotted  the  corrected  m^  terms  against  Pn  velocities. 

The  resulting  plot  is  shown  in  Figure  2.  The  line  of  Marshall  et  al  is  super¬ 
posed,  but  is  shifted  down  .12  magnitude  units  since  we  reduced  our  m^  to  that 
of  a  granite  station,  making  our  crustal  corrections  negative.  The  data  points 
are  suggestive  of  a  curve  that  is  similar  in  shape  to  that  presented  in 
Marshall  et  al.  If  averages  for  each  P^  velocity  interval  are  used,  the 
resulting  data  points  and  curve  are  not  dissimilar  to  those  in  Figure  1. 

This  shows  that  the  application  of  crustal  station  corrections  does  not  alter 
the  conclusions  of  MSR. 
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Log10(A/T)  residuals 


Am  (.(Corrected  for  crust) 


U.S 


Figure  2.  LRSM  magnitude  residuals  from  Booth  et  al  (1974)  with  crustal 

corrections  due  to  Der  et  al  (1979)  as  a  function  of  Pn  velocity 
as  from  Marshall  et  al  (1979).  The  curve  drawn  is  from  Marshall 
et  al  adjusted  down  by  0^12  to  allow  for  the  crustal  corrections. 
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Let  us  now  examine  the  worldwide  validity  of  this  relationship  by 

plotting  the  bias  teims  of  North  (1977)  against  the  values  listed  by 

MSR  (see  Figure  3).  The  bias  terms  are  not  corrected  for  the  crust,  since 

such  correction  terms  are  not  available.  The  MSR  data  set  is  dificient  in 

stations  which  have  low  P  velocity.  In  order  to  correct  this  deficiency,  we 

n 

have  added  some  African  stations  (BHA,  MTD,  KRR,  CLK  and  CIR) .  We  assume  that 
the  P^  velocity  for  these  was  7.95  as  given  by  Searle  and  Gouin  (1971).  It 
must  also  be  pointed  out  that  even  lower  P^  velocities  are  suggested  by  Dopp 
(1964)  in  the  western  rift  zone.  Again,  superposing  the  curve  of  Marshall 
et  al  resutls  in  a  strong  suggestion  that  P^  velocities  and  magnitude  biases 
are  correlated  worldwide  in  the  manner  presented  in  their  paper. 

The  U.S.  WWSSN  stations  also  show  a  correlation  between  P  velocity  and 

n 

magnitude  bias  terms  similar  to  the  LRSM  data  set  (Figure  4,  after  North,  1°77). 

In  summary,  we  can  find  no  empirical  fault  with  the  procedure  of  Marshall 
et  al.,  although  its  theoretical  foundation  is  slightly  shaky. 

Before  the  procedure  is  adapted  for  general  use  it  would  perhaps  be 
advisable  to  extend  the  existing  m^-bias,  Pn  velcoity  data  set  to  test  the 
general  validity  of  such  relationships.  In  our  view,  data  points  in  the 
foldbelts  surrounding  the  Eurasian  landmass  are  crucial  and  can  be  obtained 
without  acquiring  new  data  besides  the  WWSSN  and  the  SRO  network.  In  addition, 
other  indirect  methods  suing  other  types  of  independent  geophysical  variables, 
e.g.,  travel  time  delays  and  spectral  anomalies  also  need  to  be  tested. 

Finally,  the  defining  relationship  for  absorption  is  the  spectral  ratios  of 
common  events  at  different  stations  and  these  should  be  computed  to  the  extent 
data  are  available. 
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non  U.S. 


P„  Velocity 


Figure  3.  Magnitude  residuals  of  North  (1977)  plotted  as  a  function  of  Pn 
velocity  due  to  Marshall  et  al  (1979).  The  superimposed  line  is, 
again,  due  to  Marshall  et  al  (1979). 


pP  Corrections 

Figure  5a  from  MSR  shows  the  waveforms  resulting  from  the  superposition 
of  P  and  pP  at  different  delays.  A  more  full  set  of  superpositions  may  be 
found  in  Marshall  (1972),  from  which  presumably  Figure  5a  was  drawn.  Figure  5b 
gives  the  magnitude  corrections,  DC(T)  to  be  added  to  the  observed  magnitude 
in  order  that  the  magnitude  be  the  same  as  if  pP  had  perfectly  reinforced  the 
direct  wave.  The  correction  is  a  function  of  DCR  which  is  equal  to  the  P-pP 
delay  divided  by  the  dominant  period. 

From  a  theoretical  point  of  view.  Figure  5b  must  be  in  error  for  small 
DCR — DC(T)  must  tend  to  infinity.  The  wave  shape  will  go  to  the  derivative 
of  the  direct  pulse  in  the  limit  of  shallow  depth  (DCR  ■+■  0).  Thus  the  period 
will  tend  to  a  constant  as  DCR  -*•  0  while  the  signal  amplitude  tends  to  zero. 
Thus  the  correction  amplitude  DC(T)  should  tend  to  infinity.  As  a  practical 
matter,  of  course,  the  reflection  becomes  increasingly  imperfect  as  the 
explosion  craters  and  the  correction  must  approach  some  limit  as  is  suggested 
in  Figure  5b. 

Figure  5b  is  not  complete  in  another  respect.  It  is  produced  using  a 
single  waveform.  However,  over  the  range  of  magnitudes  and  t*'s  and  instru¬ 
ment  responses  of  interest  there  are  several  waveforms.  Figure  5c  gives  our 
calculations,  (using  techniques  discussed  in  Blandford,  1976),  of  DC(T)  for  the 
WWSSN  SP  response  and  a  granite  RDP  for  a  low  frequency  source  (Y  =  300  kt, 
t*  =  0.6)  intermediate  (Y  *  150  kt,  t*  *  0.4),  and  high  (Y  =  75  kt,  t*  =  0.2). 
The  waveforms  from  which  these  figures  are  derived  may  be  found  in  Appendix  I. 
We  see  large  DC(T)  values  for  small  DCR  as  discussed  above,  and  appreciable 
(~  ±  0.07  m^)  scatter  in  DC(T)  for  any  fixed  value  of  DCR.  Thus,  strictly 
speaking,  corrections  should  be  made  independently  for  different  instrument 
responses  and  different  "station"  t*  values. 

Further  complications  arise  when  one  notes  that  the  pP  reflection  is  not 
perfect.  Blandford  (1976)  and  earlier  workers  cited  in  that  study  showed  that 
the  typical  reflection  coefficient  in  the  range  around  1  Hz  for  NTS  explosions 
REX,  SCOTCH,  and  BENHAM  was  -0.5,  and  Shumway  and  Blandford  (1977)  showed  that 
the  reflection  coefficient  was  close  to  -1.0  at  low  frequencies.  Shumway 
and  Blandford  (1977)  also  developed  a  technique  for  producing  synthetic 


-14- 


seismograms  with  a  reflection  coefficient  which  varies  as  a  function  of 

frequency.  Using  this  technique,  we  have  repeated  the  calculations  for 

2 

Figure  5  c  but  with  a  reflection  coefficient  equal  to  -0.5(exp(-f  )  +  1)  where 
f  is  frequency  in  Hz.  The  results  are  seen  in  Figure  5d.  Here  we  see  an  even 
larger  scatter  than  in  Figure  5c.  Thus  one  concludes  that  the  MSR  correction 
is  valid  on  the  average  but  may  have  -  0.1  m^  error  in  any  particular  appli¬ 
cation. 

There  is  good  empirical  evidence  for  the  necessity  of  making  a  correction 
of  this  type.  Figures  6  and  7  from  Blandford  et  al.  (1976)  show  variations  of 
the  "a"  phase  about  the  "c"  phase  (which  is  generally  the  maximum).  Since  the  "a" 
phase  is  unaffected  by  pP  at  NTS  a  large  value  of  the  residual  (large  a)  means 
that  c  has  been  reduced  by  pP,  and  a  small  residual  (small  a)  means  that  c 
has  been  enhanced  by  pP. 

In  both  figures  the  residuals  are  plotted  versus  the  theoretical  pP  time 
(given  by  twice  the  depth  divided  by  the  uphole  average  velocity  as  determined 
from  LLL  data)  divided  by  the  dominant  period.  The  dominant  period  was 
determined  empirically  at  each  station  using  a  regression  on  period  versus 
amplitude  for  the  events  in  question.  (It  would  be  worth  doing  the  plots  over 
again  using  the  individual  period  for  each  event.)  We  see  that  effects  are 
large  and  there  is  a  great  deal  of  scatter,  in  agreement  with  the  results  in 
Figures  5c  and  5d. 

It  seems  to  us  that  a  more  soundly  based  approach  is  needed  for  this 
problem,  although  the  MSR  technique  is  better  than  nothing. 

1-  Run  the  Shumway  and  Blandford  (1977)  maximum  likelihood  estimate 
to  determine  a  and  t. 

2-  Using  these  a  and  t  values  compute  the  synthetic  waveforms  for  all 
stations — using  t*  as  appropriate — to  determine  the  m^  correction 
due  to  pP. 

Figure  8  shows  that  in  the  magnitude  range  of  interest  (-  70  -  300  kt), 
the  m^ryield  slope  is  -  0.7,  although  in  the  range  2  -  70  kt  the  slope  is 
close  to  1.0.  We  have  performed  work  to  show  that  the  analysis  techniques 
often  used  bias  observed  slopes  to  smaller  values  if  small  events  at  differ¬ 
ent  test  sites  are  included  with  large  events.  There  is  every  reason  to 
suspect  this  problem  in  the  MSR  study  at  low  yields.  The  true  slope  of  0.7 
at  large  yields  can  combine  with  the  biased  slope  <  1.0  at  lower  yields  to 
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Figure  5a. 
5b. 


(a) 

(a)  Surface  (b)  Intermediate  (c)  Deep 


(b) 


Variation  of  the  pP-P  interference  effect  with  explosion 
depth  (from  Marshall,  Springer  and  Rodean,  1979). 

The  magnitude-correction  term  DC(T)  given  as  a'  function  of 
the  depth-correction  ratio  DCR.  The  DC(T)  is  added  to 
magnitude  to  give  the  final  magnitude  m.  (from  Marshall, 
Springer  and  Rodean  (1979).  4 
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Figure  5c.  Magnitude  correction  ae  a  function  of  DCR  -  (pP  delay/c 
phase  period)  from  waveforms  In  Appendix  I.  Perfect  pP 
reflection  (qpP  ■  -1.0),  WWSSN  SP  response,  von  Seggern- 
Blandford  granite  reduced  displacement  potential.  Cor¬ 
rection  Is  to  be  added  to  achieve  maximum  mfe  possible  with 
perfect  pP  enhancement  for  each  yield. 


5d.  Magnitude  correction  as  a  function  of  DCR  ■  (pP  delay/c 
phase  period)  from  waveforms  in  Appendix  1.  Variable 
pP  reflection  (-0.5(1  +  exp(-f^)),  WWSSN  SP  response, 
von  Seggem-Blandford  granite  reduced  displacement  poten¬ 
tial.  Correction  la  to  be  added  to  achieve  maximum  m^ 
possible  with  perfect  enhancement  for  each  yield. 
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Figure  6.  RKON  station:  Residuals  of  log  (a)  with  respect  to  c  as  a 
function  of  scaled  2-way  pP  travel  time.  The  solid  line  is 
hand-drawn  as  a  subjective  estimate  of  the  main  data  trends 
(from  Blandford,  1976b). 


Figure  7.  Residuals  of  log  (a)  at  HNME  with  resoect  to  c  as  a  function 
of  scaled  2-way  pP  travel  time.  The  solid  line  is  hand-drawn 
as  a  subjective  estimate  of  the  mean  data  trends  (from 
Blandford,  1976b). 


A/T  (CORRECTED  FOR  PERIOD) 
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Figure  8.  Theoretical  amplitude-yield  curves  for  t*  -  0,  0.1,  0.2,  0.4, 

0.8;  granite,  one-half  maximum  peak-to- trough  amplitude  for  signal, 
corrected  for  LRSM  instrument  response  at  measured  period,  T, 
and  divided  by  T,  with  surface  reflection.  pP  delay  equal  to 
0.12Y1/3  sec,  with  Y  in  kt.  From  Blandford  (1076b). 
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give  an  (incorrect)  smooth  straight  line  of  slope  ~0.7  from  low  to  high 
yields.  Note  that  there  are  discontinuities  in  the  m^:yield  curve  around 
100  kt  due  to  the  effects  of  pP  in  introducing  "false  cycles"  of  varying 
size.  This  suggests  that  in  determining  m^,  at  shield  sites  especially, 
the  period  correction  should  not  be  made.  (This  figure  is  for  the  LRSM 
response.  For  the  WWSSN  response  which  peaks  at  lower  frequencies,  the 
"false  cycles"  are  not  as  great  a  problem.) 

Figure  9  plotted  from  waveform  values  in  Appendix  I  shows  directly 

the  variation  of  log(c/GT)  as  a  function  of  log  yield  for  the  granite 

2 

reduced  displacement  potential,  t*  =  0.4,  apP(f)  =  -0.5(1  +  exp(-f  )).  We 
see  that  several  magnitude : yield  slopes  are  possible,  depending  on  the 
variation  of  explosion  depth  (pP  delay)  with  yield.  Constant  depth  gives 
a  slope  in  the  range  0.7-0. 8.  If  depth  varies  as  the  cube  root  of  the 
yield  the  pP  delay  would  vary  from  0.1  sec  at  75  kt  to  0.16  sec  at  300  kt 
resulting  in  a  slope  of  0.83.  Slopes  of  1.0  or  greater  are  possible  with 
different  depth  relations.  Assuming  that  the  pP  >_  0.1  sec  then  the  possible 
range  in  magnitude  is  about  0.15-0.2.  For  perfect  pP  reflections  as  seen  in 
Figure  10  the  range  is  larger  0.25-0.3,  but  there  is  no  good  evidence  of 
which  we  are  aware  for  such  high  reflection  coefficients  at  1  Hz. 

As  can  be  seen  in  Figure  11,  the  variation  of  m^  due  to  plausible 
variation  in  t*  is  much  larger,  approximately  0.6  magnitude  units.  The 
variations  of  different  measures  of  magnitude  are  in  parallel  as  a  function 
of  t*  and  are  substantially  unchanged  whether  measurements  are  made  on  a,  c, 
c  corrected  for  period :cg,  or  on  c  corrected  for  period  and  divided  by 
period :m£  •  log(c/GT). 

In  Figure  11  we  can  see  that  there  is,  accompanying  this  variation,  a 
variation  which  is  basically  undetectable  on  WWSSN  film  of  the  c  phase  period 
from  0.8  to  0.9  seconds.  The  a  phase  period  varies  more,  from  0.5  to 
0.8  sec;  however  the  period  of  the  first  half-cycle  is  very  hard  to  measure,  and 
the  particular  parameter  choice  selected  to  be  plotted  has  overemphasized  the 
variation,  which  is  more  typically  0.6  to  0.8  sec.  From  this  figure  one  may 

correction  of  6m^  *  1.5  6t*  which  nicely  falls  near  the  "spectral" 
estimates  Sm^  ■  (ir/loge10)f 6t*  which  gives  Sm^  *  1.376t*  for  T  *  1,  and 
-  1.705t*  for  T  ■  .8. 

For  convenience,  we  have  included  Figures  12  and  13  which  give  the 
amplitude tyield  relations  for  the  a  and  c  phases  for  perfect  reflection. 
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log  Y 


Figure  9.  Theoretical  amplitude:yield  for  t*  =  0.4,  granite;  c  phase 

corrected  for  WWSSN  response  at  measured  period  T,  and  divided 
by  T,  with  surface  reflection  coefficient  apP(f)  ■  -0.5(1  + 
exp(-f2)),  pP  delay  equal  to  .05,  .1,  .15,  .2,  .3,  .4,  .ft  sec 
for  Y  ■  75,  150,  300  kt.  From  waveforms  in  Appendix  I. 
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Figure  in.  Theoretical  amplitude  .’yield  for  t*  ■  0.4  granite;  c  phase  corrected 
for  WWSSN  response  at  measured  period  T,  and  divided  by  T,  with 
surface  reflection  coefficient  apP(f)  ■  -1.0  pP  delav  equal  to 
.05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y  *  75,  150,  300  kt.  From 
waveforms  in  Appendix  I. 


DISCUSSION  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


It  seems  that  corrections  for  effects  of  pP  and  absorption  should  be 
made  in  yield  estimation  and  that  the  procedures  outlined  in  MSR  are  rea¬ 
sonable  approaches  to  the  problem.  The  pP  problem  could  be  attacked  in  a 
more  operational  manner  and  with  more  accuracy  by  determining  the  pP  delay 
and  reflection  coefficient  and  then  computing  the  correction  directly  from 
the  appropriate  synthetic  waveform.  The  absorption  problem  should  be 

attacked  by  further  analysis  of  the  P  -residual  correction. 

n 

Another  approach  to  the  absorption  problem  which  seems  very  promising  Is 
that  of  determining  t*  from  the  periods  of  short-period  S  waves  of  deen  events 
with  plenty  of  high-frequency  energy.  This  approach,  in  contrast  to  the  mag¬ 
nitude  of  travel-time  residual  methods,  does  not  require  that  the  precise 
event  location  or  magnitude  be  known  or  even  that  the  absolute  calibration 
level  of  the  instrument  be  known.  WWSSN  data  and  even  data  reported  from  the 
AI,  can  be  used  in  these  determinations.  The  t*  is  determined  in  this  approach 
exactly  in  the  1  Hz  frequency  band  that  is  of  greatest  interest  for  magnitude 
determinations . 

Other  sources  of  error  in  magnitude  determinations  are  the  focusing- 
defocusing  scatter  in  magnitude  generated  by  small  (20-50  km)  shifts  in  the 
test  site,  effects  of  clipping  and  of  signals  below  the  noise  level,  ampli¬ 
fication  due  to  crustal  structure  at  the  recording  site  and  at  the  source  and, 
of  course,  variations  in  coupling  due  to  the  shot  point  medium. 

It  could  be  a  worthwhile  study  to  determine  magnitudes  for  selected 
shots  of  interest  including  PILEDRIVER,  HARDHAT,  RUBIS,  SAPHIRE,  SALMON, 
and  selected  Soviet  shots  with  all  of  the  above  considerations  taken  into 
account.  The  principal  data  base  should  be  the  WWSSN  system  with  LRSM 
supplementation  for  the  HARDHAT-PILEDRIVER  pair,  because  HARDHAT  occurred 
before  much  of  the  WWSSN  system  was  up.  AEDS  data  should  also  be  analyzed 
in  order  to  estimate  and  allow  for  network  biases  and  inescapable  scatter. 
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Figure  11.  Theoretical  curves  of  magnitude  and  period  as  a  function  of  t*  for 
several  definitions  of  magnitude.  Granite,  WWSSN  response  surface 
reflection  coefficients  apP(f)  »  -0.5(1  +  exp(-f^)),  pP  delay 
equals  0.15  sec,  Y  -  150  kt.  Cg  is  c  corrected  for  instrument 
response  at  period  T,  c/G,  and  mc  ub  c/GT.  From  waveforms  in 
Appendix  I. 
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Figure  12.  Theoretical  amplitude ’.yield  curves  for  t*  *  0.4  granite,  a.  phase 
with  surface  reflection  coefficient  apP(f)  -  -1.0,  pP  delay  equal 
to  .05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y  -  75,  150,  300  kt. 

From  waveforms  in  Appendix  I. 
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Figure  13.  Theoretical  amplitude:yield  curves  for  t*  ■  0.4  granite,  £  phase 
with  surface  reflection  coefficient  opP(f)  ■  -0.1,  pP  delay  equal 
to  .05,  .1,  .15,  .2,  .3,  .4,  .6  sec  for  Y  -  75,  150,  300  kt.  From 
waveforms  in  Appendix  I. 
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APPENDIX  I 

WWSSN  Waveforms  and  Short-Period  Measurements  of 
75,  150,  300  kt  granite  explosions,  t*  *  0.2,  0.4, 
0.6,  pP  delay  .05,  0.1,  0.15,  0.2,  0.3,  0.4,  0.5,  0.6 
seconds.  Full  and  partial  pP  reflection  coefficients. 


A-l 


Appendix  I 


This  appendix  contains  theoretically  calculated  waveforms  together  with 

measurements  made  from  those  waveforms  by  computer.  These  waveforms  serve  as 

the  basis  for  many  of  the  plots  in  the  body  of  this  report.  The  calculation 

technique  is  to  form  an  amplitude  spectrum  as  the  product  of  the  von  Seggern- 

Blandford  scaled  reduced  displacement  spectrum  for  granite  times  the  WWSSN 

short-period  instrument  response  and  finally  times  exp(-mt*/2) .  The  minimum 

phase  wavelet  with  this  spectrum  is  then  convolved  with  two  delta  functions 

representing  P  and  pP.  In  the  case  where  pP  has  the  spectrum  modified  by 
2 

-0.5(1  +  exp(-f  )),  the  second  delta  function  is  replaced  before  convolution 
by  the  minimum  phase  wavelet  with  that  spectrum. 

The  first  set  of  waveforms  is  for  perfect  pP  reflection  (apP  =  -1.0) 
and  the  second  set  apP(f)  =  -0.5  (1  +  exp(-f^)).  Yields  of  75,  150,  and  300 
kt  are  put  through  t*  values  of  0.2,  0.4,  and  0.6.  Both  the  a  and  c  phases 
are  measured.  For  each  waveform  the  upper  set  of  four  numbers  is  period  (T) , 
log  amplitude  (logA) ,  log(A/G)  where  G  is  the  gain  relative  to  1  Hz  at  period 
T,  and  logA/(GT)  which  is  proportional  to  magnitude  as  usually  defined. 

The  upper  set  of  four  numbers  is  for  the  a  and  phase  the  lower  set  is 
for  the  c  phase.  All  the  magnitudes  are  correct  relative  to  each  other; 
however,  the  overall  absolute  level  has  been  adjusted  to  be  in  the  range 
of  typical  values. 
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